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We study structure of excited states of ^"Be with the method of variation after spin parity 
projection in the frame work of antisymmetrized molecular dynamics. Present calculations describe 
' many excited states and reproduce the experimental data of E2 and El transitions and the new 

, data of the /3 transition strength successfully. We make systematic discussions on the molecule-like 

^\ • structures of light unstable nuclei and the important role of the valence neutrons based on the results 

T-H ' obtained with the framework which is free from such model assumptions as the existence of inert 

cores and clusters. 
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I ■ I. INTRODUCTION 

(N : 

I I Owing to the radioactive nuclear beams, the data of unstable nuclei are increasing rapidly. Recently the structure 
• of excited states as well as the ground states are very attractive in the study of unstable nuclei. In the light nuclear 
00 ' region, one of the hot and important subjects has been the clustering features. The clustering structures of the ground 
, states are already well known to be developed in some light ordinary nuclei such as ^Li, ^Be and ^°Ne. In the studies 
of the light unstable nuclei the clustering structures were predicted also in the very neutron-rich nuclei. 

1^ , It is natural that various molecule-like states may appear in the excited states of light unstable nuclei because the 
excitation due to relative motion between clusters is important in the light nuclear region. The clustering structure 
must be one of the important aspects to understand the exotic features of the unstable nuclei. However there remain 
many mysterious problems of clustering development and the role of valence nucleons in the unstable nuclei. Our 
I aim is to make systematic study of structure change with the increase of the excitation energy. We are to find new 
I , . features of unstable nuclei such as exotic shapes and molecule-like structures in the excited states and to understand 
^ the mechanism of the various structures. 

^ ^°Be, one of the challenges in the study of light unstable nuclei, were and are studied hard theoretically |§,^||5]j(],0.|[ . 

Recent experiments of the charge exchange reactions ^°B('^He,t)^'^Be [|j let us know the strength of the Gamow- Teller 
J> [ transitions to the excited states of ^°Be. These new data of f3 transition strength which are deduced from the cross 
k>( " sections at the forward angle are very helpful to study the structure of the excited states. 

rN ] In the theoretical researches on the structure of the ground states of Be isotopes such as the reference [|| , it was 
. found that the developed 2a clustering structure in *Be weakens in the ^'^Be because of the excess two neutrons. As 

" " ' for the excited states of ^"Be, the previous calculations Jl3] with a simplest version (variation after parity projection) 
of antisymmetrized molecular dynamics (AMD) Pp pc| , p^ predicted that the largely deformed states with the 2Q;-|-2n 
clustering structure construct a rotational band K — 1^ in low energy region which is supported by the experimental 
energy levels. 

The problem of the abnormal spin parity 1/2+ in the ground state of ^^Be is considered to have a close relation 
with the deformed non-normal parity states of ^°Be. In the simple shell model consideration the normal parity of ^^Be 
is negative because of the valence neutron in the pi/2 shell, however, the ground states is known to be J'^ = 1/2+. 
It is suggested that one of the reasons for the parity inversion in ^^Be is the energy gain of the sdi/2 mixing orbit 
in the prolately deformed system due to the developed clustering. The s(ii/2 orbit of neutrons in neutron-rich Be 
isotopes is interesting and essential to understand the structure of light neutron-rich nuclei and is easily connected 
with the deformed structure of the non-normal parity states of ^"^Be. Furthermore it is natural that the negative 
parity band due to the deformed states of ^"Be makes us imagine a corresponding positive parity band constructed 
by the analogous deformed intrinsic state where the neutrons in the sdi/2 orbit play important roles. 

Recently W. Von Oertzen proposed a dimer model |^ to describe the excited states of Be and B isotopes system- 
atically. In his idea he supposed 2a clusters and surrounding nucleons in molecular orbits like a and tt bonds as the 
orbits of electrons in the molecule. The neutron-rich Be nuclei can be written by 2a and neutrons which occupy the 
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molecule-like a and tt orbits. This idea is helpful to understand the level structure of many excited states of ^Be and 
^"^Be and to study the mechanism of the deformation of Be isotopes. 

A. Dote et al. studied Be isotopes by microscopic calculations with the simplest version of AMD under a constraint 
on the deformation parameter |5[|. They analyzed the single-particle orbits of the valence neutrons and found the 
importance of molecular orbits. It was found that in the negative parity state of ^''Be the main component of the 
valence neutrons is the positive parity orbit which seems to be the a orbit in terms of the dimer model. They tried 
to represent the excited state of ^"Be by constructing 2p-2h states using the wave function of the negative parity 
state which is considered to be Ip-lh state. Comparing the optimum deformation parameters for three states 0^, 1~ 
and Oj, they found that the intrinsic deformation is larger in 1^ than 0^ and largest in Oj. The results is consistent 
with the ones predicted by W. Von Oertzen. They calculated excited states 0^ of ^"Be without assuming any cluster 
cores for the first time and confirmed that the 2a cluster cores may exist in the excited states. However since the 
total-angular-momentum projection was not done microscopically, the discussion didn't go beyond the intrinsic system 
and it is not easy to study other excited states. They have not mentioned the data such as E'l, E2 and (3 transitions 
which are important in the study of excited states. 

^°Be has been studied with other microscopic models such as multi cluster models 0, AMD+GCM(Generator 
Coordinate Method) joj and the parity-projected Hartree-Fock(HF) calculations In the former two method the 
existence of a clusters was assumed. In the last work with parity-projected HF calculations they assumed the excited 
0^ state as a 2p-2h state by choosing the positive parity eigen state for the valence two neutrons. In the work with 
HF calculations, the analysis was limited in the intrinsic system and it is not easy to study the higher excited states. 

Our aim is to make systematic studies on the structure of many excited states of the light unstable nuclei with a 
microscopic model which is free from such assumptions as the inert cores, the existence of clusters, and the particle- 
hole configurations. For this aim we adopt a method of antisymmetrized molecular dynamics (AMD VAMD has been 
already proved to be a very useful theoretical approach for the structure of the light nuclei P,pp| JTl|]T^ . In the AMD 
framework basis wave functions of the system are written by Slater determinants where the spatial part of each single- 
particle wave function is a Gaussian wave packet. In the previous works on the structure of Li, Be and B isotopes , 
we applied the simplest version of AMD in which the energy variation is done after the parity projection but before 
the total-angular momentum projection. It was found that the structure of the ground state changes rapidly between 
the shell-model-like structure and the clustering structure as the neutron number increases up to the neutron-drip 
line. The AMD calculations succeeded to reproduce well the experimental data for electro-magnetic properties such 
as the magnetic dipole moments and the electric quadrupole moments. It owes to the flexibility of the AMD wave 
function which can represent the clustering structures of light unstable nuclei without assuming any inert cores and 
clusters. 

In the present paper, we study the structure of the excited states of -'^^Be by performing variational calculations 
after the spin parity projection with finite-range interactions in the framework of AMD. For higher states, AMD wave 
functions are superposed so as to be orthogonal to the lower states. The AMD approach of the variation after spin 
parity projection(VAP) has been already found to be advantageous for the study of excited states of the light nuclei 
and be very useful to describe various molecule-like structures [|l3|,0 . We do not rely on any assumptions such as the 
inert cores, the single-particle orbits in the mean field, and the existence of clusters. By microscopic calculations of the 
expectation values of the corresponding operators, we can acquire the theoretical values of £'2, El and (3 transitions 
which are the valuable informations of the excited states to be compared directly with the experimental data. 

In Sec. H, we explain the formulation of AMD for the study of the nuclear structure of excited states. The effective 



interactions are described in Sec. Ill, and the results are presented in Sec. IV comparing with the experimental data. 
In Sec. 0, the detailed discussions on the structures and on the mechanism of the development of the molecule-like 



structures are made. Finally summary is given in Sec. VI 



II. FORMULATION 

In this section we explain the formulation of AMD for the study of nuclear structure of the excited states p^Jl^ . 

A. Wave function 

The wave function of a system is written by AMD wavefunctions, 

^^c^amd+c^'amd + ---- (1) 
An AMD wavefunction of a nucleus with mass number A is a Slater determinant of Gaussian wave packets; 
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'^AMoi'Z) = ^A{ipi,(p2,---,VA}, (2) 



where the i-th single-particle wave function ipi is a product of the spatial wave function, the intrinsic spin function 
and the isospin function. The spatial part (jjXi is presented by complex parameters Xn, X2i, X^i, x^i is the intrinsic 
spin function parameterized by and is the isospin function which is fixed to be up(proton) or down(neutron) 
in the present calculations. Thus an AMD wave function is parameterized by a set of complex parameters Z = 
{Xni,(,i} {n — 1,3 and i = 1,^), where X^'s indicate the centers of Gaussians of the spatial part and ^i's are the 
parameters for the directions of the intrinsic spins. 

If we consider a parity eigen state projected from a AMD wave function the total wave function consists of two 
Slater determinants, 

*(Z) = (l±P)$AAfi)(Z), (4) 

where P is a parity projection operator. In case of a total-angular momentum eigen state the wave function of a 
system is represented by integral of the rotated states, 

$(Z) = P-ltK''^AMDm = J dnDi^K'{^)R{^)'^AMDiZ), (5) 

for which the expectation values of operators are numerically calculated by a summation of mesh points on the Euler 
angles fl. 

In principal the total wave function can be superposition of independent AMD wave functions. In order to construct 
higher excited states we consider superposition of the spin parity projected AMD wave functions P^j%,^AMD, 

$ = cP/,± ,$amd(Z) + c'P^,%<^>amd{Z') + ■■■. (6) 



The detail is mentioned later in subsection [I E 



B. Energy variation 

We make variational calculations to find the state which minimizes the energy of the system; 

{mm 
($1$) 



(7) 



by the method of frictional cooling. Regarding the frictional cooling method in AMD, the reader is referred to papers 
For the wave function ^(Z) parameterized by Z, the time development of the parameters are given by the 
frictional cooling equations, 

d^k , . , 1 d (ci>(Z)|g|<i>(Z)) 

with arbitrary real numbers A and /i < 0. It is easily proved that the energy of the system decreases with time. After 
sufficient time steps of cooling, the parameters for the wave function of minimum-energy state are obtained. 
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C. Lowest J states 



In order to obtain the wave function for the lowest J state, we perform the energy variation for the spin parity 
eigenstates projected from an AMD wave function. In this case the trial function is $ = P//^,$yiMD(Z). In the 
previous works with the simplest version of AMD for the study of nuclear structure, the approach was variation 
after only the parity projection but variation before the total spin projection(VBP). In the present paper, the approach 
is variation after the spin parity projection(VAP). First we make VBP calculations to prepare an initial wave function 
^AMDC^init) for the VAP calculations. We choose an appropriate K' quantum that gives the minimum diagonal 
energy of the spin parity eigenstate (Pj;^'^,^'AMi)(Z„it)|i?|P//^,$AMD(Zi„it))/(P//^,(Z„it)|P;^^,(Zi™t)), where 
K' is the component of the total-angular momentum along the approximately principal axis in the intrinsic system. 
For each spin parity J±, we perform VAP calculation for {P^j\,<^>amd{Z)\H\p1^j^,^amd{Z)) / {P^%{Z)\P^^j^,{Z)) 
with the adopted K' quantum from the initial state. In the VAP procedure, the principal z-axis of the intrinsic 
deformation is not assumed to equal with the 3-axis of Eulcr angle in the total-angular momentum projection. In 
general the principal z-axis is automatically determined in the energy variation, that is to say that the obtained 
state P^f^i^AMD by VAP with a given K' = (J3) can be the state with so-called if- mixing in terms of the intrinsic 
deformation where the K = (J^) quantum is defined on the principal axis. In many cases, z-axis determined in the 
VAP calculation is found to be approximate to the 3-axis. It means that the obtained states do not contain if -mixing 
so much, and K' is considered to be the approximate K quantum for the principal axis. The excited states in the 
band = K"^ other than K' quantum of the lower J* states are obtained by VAP calculations for P^.^,,^amd 
under the constraint on the principal axis z as equal to the 3-axis. 



D. Higher excited states 

As mentioned above, with the VAP calculation for ^(Z) = P//^/$amu(Z) of the eigen state with K' , we obtain 
the set of parameters Z = Zj'* which present the wave function for the first J* state. To search the parameters Z 
for the higher excited J* states in the n-th = K'^ band, the wave functions are superposed to be orthogonal 
to the lower states as follows. The parameters Z;[^ for the n-th J* state are reached by varying the energy of the 
orthogonal component to the lower states; 

if,(7\ pJ± ffy (7\ {PAm'^AMD{Zi^)\P-li\,<i>AMD{Z)) j± j± 

fc-1 {Pmk''^amd[Z^ }\Pj^j^,^amd[Z^ )} 



E. Expectation values 

After VAP calculations for various states, the intrinsic states ^amd^ ^\mD' ' ' ' ' ^amD' which approximately 
correspond to the states, are obtained as much as the number of the calculated levels. Finally we determine the 
wave functions for the states by diagonalizing the Hamiltonian matrix {P^j^' ^am IPaik" ^am d) ^^'^ norm 
matrix {Pajk'^amd\Pmk"^amd) simultaneously with regard to for all the intrinsic states and {K' ,K"). In 
comparison with the experimental data such as energy levels and E2 transitions, the theoretical values are calculated 
with the final states after diagonalization. 



III. INTERACTIONS 



The adopted interaction for the central force is the case 3 of MVl force jl^] , which contains a zero-range three-body 
force V'^^^ as a density dependent term in addition to the two-body interaction V'^'^^ of the modified Volkov No.l force. 



Vdd = V^^^ + 



= (1 



hP^ - hPr - mP^Pr) { Va cxp 



Vfjexp 



if-:)- 



Va = -83.34MeV, = 1.60fm, Vr = 99.86MeV, = 0.82fm, 
^(3) ^ ^;(3)j(i.^ _ Y2)5{vi - ra), = SOOOMeV fm^, 



(11) 
(12) 

(13) 
(14) 



4 



where Per and Pt stand for the spin and isospin exchange operators, respectively. As for the two-body spin-orbit force 
Vls, we adopted the GSRS force Jlq] as foUows, 



Vls = {^^/exp {-Kir^)+uiie,,p (-Kur^)} Il±^Il±^L • (Si + S2), 



Ki = S.Ofm"^ Kii = 2.778fm" 



(15) 
(16) 



The Coulomb interaction Vc is approximated by a sum of seven Gaussians. 

We treat the resonance states within a bound state approximation by situating an artificial barrier £barr- In the 
variational calculation we variate the modified energy £' = £ + f barr with the additional artificial energy f barr instead 
of = (<I>|i7|<I>)/(<I>|$), where H is the Hamiltonian operator; H = T+Vdd + Vls + Vc- For the barrier in the present 
calculations we adopt a function of the real part of the centers of the single-particle wave functions as follows, 



£ha 



Vw X] *^^P 



-21/ r. 



\Rc\Xi 



(17) 



=ro(2Max[Z,iV])'/V6 (18) 
Vw = 2.5 MeV, ro = 1-2 fm, b = ^ fm, (19) 

where Z and N are proton and neutron numbers, respectively. We put the artificial barrier only in the variational 
calculation, but do not put the barrier in the diagonalization of the Hamiltonian matrix after VAP calculations. In 
the present calculations of ^°Be, the artificial energies due to the barrier are found to be a few hundreds keV at most 
and make no significant effect on the results in most levels. 



IV. RESULTS 



The structure of the excited states of -'^^Be is studied with the VAP calculations in the framework of AMD. In 
this section we display the theoretical results of the excitation energies, E2, El, and /3 transitions which can be 
directly compared with the experimental data. The detail of the structures are discussed in the next section. The 
adopted parameters of the interactions are m = 0.62, b = h = ior the Majorana, Bertlett and Heisenberg terms of 
the central force and the strength of the spin-orbit force uj = —un — 3000 MeV (case(l)). Trying another set of 
parameters case(2) with m = 0.65, b = h = and uj — ~ujj — 3700 MeV, we did not find significant differences 
in the results. The set of parameters of case(l) is the one adopted in the work on ^^C On the other hands, 

the VAP calculations with the set of interactions case(2) reproduce the parity inversion of the ground state of ^^Be. 
The optimum width parameters v of wave packets are chosen to be 0.17 fm^^ for case(l) and 0.19 fm for case(2) 
which give the minimum energies in VBP calculations of ^"^Be. The resonance states are treated within a bound state 
approximation by situating an artificial barrier out the surface as mentioned in section pH. 

The lowest J± states are obtained by VAP calculations for /'//^,$aa/d with {J'^,K^ (0+,0), (2+,0), (3+,+2), 
(4+,0), (1",-1), (2~,-l), (3~,-l), (4~,-l). Considering Oj state to be 0+ state in the second K'' = 0+ band, 
the O2 sta te is calculated by VAP as the higher excited state orthogonal to the lowest 0^ state as explained in the 
subsection ^ID| . That is to say that the Oj state is obtained by VAP for $(Z) in Eq. |l|with {J^,K', n) = (0+, 0, 2). In 
the case of higher 2"*" states, we put the constraint that the approximately principal z-axis of the intrinsic deformation 
equals to the 3-axis of the Euler angle in the total-spin projection. According to VBP calculations the second 2+ state 
is described as the band head of the lowest = 2+ band. Therefore we construct the state 2^ by choosing ( J^, K') 
of Plj^i^AMD as {J^,K') = (2'*", -1-2) under the constraint on the principal z-axis which keeps the approximately 
orthogonality to the lowest 2+ state with (J*, K') — (2+, 0). The third 2^ state is easily conjectured to be 2+ state in 
the second K"" = band like Oj state. We obtain the 2^ state by VAP for $(Z) in Eq.^with {J^,K', n) = (2+, 0, 2) 
in the same way as the Oj state, under the constraint on the principal z-axis mentioned above. It mean s that the 
orthogonal condition to 2^ is kept by superposing two wave functions as described in the subsection [ID , while the 
orthogonality to 2^ (K'^ = 2^) is taken into account by choosing the different K quantum K' ^ 0^ . 

The binding energy obtained with the case(l) interactions is 61.1 MeV, and the one with case(2) is 61.3 MeV. The 
excitation energies of the results are displayed in Fig.^. By diagonalization of the Hamiltonian matrix the excited 

lal band K'^ = OJ and 5~ is seen in the = 1~ band. Comparir 



states 4^, 6+ are found in the rotational band K'^ = Oj and 5~ is seen in the = 1~ band. Comparing with the 
experimental data, the level structure is well reproduced by theory. Although it is difhcult to estimate the width of 
resonance within the present framework, the theoretical results suggest the existence of 3+, 4+, 6+ and 5^ states 
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FIG. 1. Excitation energies of the levels in ^°Be. Theoretical results are calculated by the diagonalization of the states 
obtained with VAP by using the interactions case(l). 
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TABLE I. E2 and El transition strength. The theoretical results with the interactions case(l) are compared with the 
experimental data [[15[. 



transitions 


Muh. 


exp. 


theory (e fm'^) 


2+^0+ 


E2 


10.5±1.1 (e fm^) 


11 (e fm^) 


0+ ^ 2+ 


E2 


3.3±2.0 (e fm^) 
1.3±0.6xl0"^ (e fm) 


0.6 (e fm^) 
0.6x10"^ (c fm) 


0+ ^ 1- 


El 



which are not experimentally identified yet. The excited levels can be classified as the rotational bands K'^ = OJ", 2+, 
and 1^. The intrinsic structures of these bands are discussed in detail in the next section. 

The data of transition strength are of great help to investigate the structures of the excited states. The theoretical 
results with the interaction case (1) and the experimental data of E2 and El transition strength are presented in 
Table |. The theoretical values well agree with the experimental data. 

The strength of the (3 decays of Gamow-Teller(GT) type transitions can be deduced from the cross sections at 
the 0° forward angle of the charge exchange reactions which have been measured recently These new data 
for the Gamow- Teller type (3 transitions are very useful probes to discuss the structures of the excited states of 
unstable nuclei. Table || shows the values of B{GT). The experimental values for the (3 transitions from ^°B(3+) to 
^"Be* are deduced from the data of the reaction ^*'B(t,'^He)^°Be. As for the theoretical values, the wave functions 
for the neighbor nucleus ^°B are calculated with VAP where (J*, if') are chosen to be (S"*",— 3) for the ground 
3^ state and (1+,— 1) for the 1^ state. ^°Be and ^°B are calculated with the case (1) and (2) interactions. The 
theoretical values reasonably match to the experimental data. The strength is not so sensitive to the interactions 
except for the decay i°B(3+) ^i°Be(2+). The results of the GT transition i°B(3+)-^i*'Be(2+) with case(l) and 
case(2) interactions underestimate the experimental data, while with another set of interactions; Volokov No.2(m = 

0.38,6 = 0,2,/i = -0.4)+G3RS(w/ = -u// = 1600MeV)+Cloumb B{GT) is calculated to be 0.41 (f|) which is 
larger than the experimental data. Since the data for ^"6(3+) ^^*'Be(9.4McV) well correspond to the theoretical 
value of -'^''B(3"'") — >^°Be(3f ), it is natural to consider the excited level of ^'^Be at 9.4MeV as the 3^ state. The result 
of B{GT) for -'^°Be(0^) ^^°B(1+) is consistent with the experimental data of the /? decay from the mirror nucleus 
i°C(0+) ^i°B(l+). It is suggested that the B{GT) for i°B(3+) ^^°Be{2j) is not large because of the well developed 
clustering structure of ^*'Be(2;^). We will discuss again the (3 decays combining with the structures of excited states 
in the next section. In the reference Q of the measurements of the strength B(GT), they mentioned about the isospin 
symmetry violation of /3 decay strength comparing B(GT_) for ^°B(3"'") ^^°Be(2j) with B(GT+) for the mirror 
reaction to ^°C. In the present calculations, we could not find the difference between B(GT_) and B(GT+) because 
there is little violation of isospin symmetry between the structures of ^°Be(2j) and the mirror state of ^"C. Taking 
into account of the experimental imformation that there exists a significant difference between the excitation energies 
of the daughter states of ^"Be and ^"C, the isospin violation of the structures of these mirror states may be possible. 
We need more detailed investigations to discuss the problem of the isospin symmetry violation of B(GT). 

V. DISCUSSION 

In this section we discuss the structure of the excited states by analyzing the wave functions. Even though the states 
obtained by VAP mix after the diagonalization of the Hamiltonian matrix, the state P//^,<I>/im_d(Z;{*) projected from 
a Slater determinant obtained in VAP with ( J^, K' , n) is dominant in the final result oi^ the Jjf state. In this section 
we consider the Slater determinant $^md(Z;^^) as the intrinsic state for the state. 

A. intrinsic structure 

In the excited states, various kinds of structures are found. Here we analyze the structures of the intrinsic states 
'^AMoi'^n'^)- It is found that the excited levels are classified into rotational bands as Of, 2f , 4f states in K'^ = 
band, 2j, 3^" in K'^ = 2+ band, Oj, 2'^, 4^, 6+ in the second K"" = Of band and 1", 2", 3", 4", 5" in K'" = 1" 
band. Particularly the molecule-like states with the well-developed 2a cores construct the rotational bands — Of 
and 1~ in which the level spacing is small because of the large moments of inertia. The density distributions of matter, 
protons and neutrons in the intrinsic states $^mli(Z;J*) are presented in Fig. 0. We found the 2a + 2n structures in 
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TABLE II. B{GT) of /3 decays with the unit The experimental data are the values^"-* deduced from the cross sections 
of ^°B(t,^He)i"Be* at 0° forward angle § and the one from the logft value jl6|. The theoretical values are obtained with 
case (1) and case (2) interactions for ^'^Be and ^'^B. 





cxp. 




initial 


final 




{r,E^) (MeV) 


{r,E^) (MeV) 


B(GT) (ll) 


^"B(3+,0) 


'"Be(2+,3.37) 


0.13± 0.05"' 


ioB(3+,0) 


i«Be(2+,5.96) 


1.51± 0.21"' 


i°B(3+,0) 


i°Be{2+ or 3+, 9.4) 


0.49± 0.13"' 


i''C(0+,0) 


i°B(l+, 0.72) 


6.25'"' 


theory case(l) 


"'B(3+) 


^"Be(2+) 


0.02 




'°Be(2+) 


1.7 




iOBe(3+) 


0.61 




iOBe(4+) 


0.12 




i«Be(2+) 


0.05 


i°Be(0+) 


"'B(l+) 


4.4 


theory case(2) 


^"B(3+) 


^"Be(2+) 


0.00 




'°Be(2+) 


1.4 




iOBe(3+) 


0.58 




iOBe(4+) 


0.14 




'°Be(2+) 


0.01 


"Be(0+) 


ioB(1+) 


3.7 



most of the intrinsic states. The density of protons indicates that the clustering structure develops more largely in 1^ 
than 0]^ and most remarkably in 0^. As seen in Fig. ||the intrinsic structure of state has an axial symmetric linear 
shape with the largest deformation, while 1~ state has an axial asymmetric shape because of the valence neutrons. 
The structures of 0^ and Ij" states are similar to the ones of the previous work with the simplest version of AMD 
||12| . The increase of the degree of the deformation along Q'l , 1^ and 0^ is consistent with the previous works such 
asP,^|j. In the ~ \~ band, the deformation toward the prolate shape shrinks as the total spin J increases. 

In the = band, the 2a cores weaken with the increase of the total spin due to the spin-orbit force. The 
reduction of the clustering structure is more rapid in the case of interaction (2) with the stronger spin-orbit force and 
the 2]*" and A'^ states in the case (2) interactions contain the dissociation of a. Regarding the dissociation of the a 
cores, the structures of those states are sensitive to the strength of the spin-orbit force. 

The 2^ has the remarkably developed molecule-like structure with an exotic shape and belongs to the rotational 
band if = . On the other hand, the lowest 3+ state of is the ordinary state in terms of p-shell. The 
theory predicts a small value of the Gamow- Teller strength ^°B(3+) ~>^''Be(23 ) because of the difference between the 
structures of the initial and final states. 



B. behavior of valence neutrons 

Even though we did not assume the existence of any clusters in the model, we have found the 2a + 2n structures 
in most of the intrinsic states as mentioned above. We study the behavior of the valence neutrons surrounding 2a 
by analyzing the single-particle wave functions to understand the role of the valence neutrons in the neutron-rich 
Be nuclei. Considering that an intrinsic state is written by a Slater determinant ^amd, the single-particle wave 
functions and the single-particle energies of an intrinsic state are determined by diagonalizing the single-particle 
Hamiltonian by the analogy with Hartree-Fock theory. First we transform the set of single-particle wave functions ipi 
of an AMD wavefunction into an orthonormal base (pa ■ The single-particle Hamiltonian can be constructed by use of 
the orthonormal base as follows ; 
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FIG. 2. The density distribution of matter, protons and neutrons of the intrinsic states are shown at left, middle and right, 
respectively. The density is integrated along the axis perpendicular to adequate planes. The figures are for the results with the 
interactions case(l). 
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FIG. 3. Single-particle energies in the intrinsic system of the 0^, 1 and states. The energies of protons(neutrons) are 
displayed in the left (right) side in each figure. 

A 

haf) = {(pa\i\(pi}) + ^((^qI^7 I f) 1(^/5(^7 - 'f'-f'f'p) 
7 

1 ^ 

7,(5 

where the Hamiltonian operator is written by a sum of kinetic term, two-body interaction term and three-body 
interaction term; H ^J2ii + J2i<j *2 + J2i<j<k ^3- 

The single-particle energies in the Oj^, 1^ and 0^ states are shown in Fig. ^. In each state the neutrons in the four 
orbits from the bottom correspond to the neutrons in the 2a clusters. The level spacing of these four lower orbits 
becomes smaller in 1^ than in 0^ and smallest in the Oj state with the increase of the distance between clusters. We 
consider the last two neutrons in the higher orbits as valence neutrons surrounding 2a cores. We display the density 
distributions of single-particle wave functions for the two valence neutrons in the left column of Fig. ^. Figures in 
the middle and right columns of Fig. ^ are for the normalized density of the positive and the negative parity eigen 
states projected from the single-particle wave functions, respectively. By analyzing the single-particle wave functions 
it is found that two valence neutrons of the Oi states contain the negative parity components more than 80% (the 
right column in figure ^ which seem to be n bonds(Fig. ^(a)) in terms of molecular orbits. On the other hand in 
the case of 0^ band, the last 2 neutrons are predominantly in the positive parity orbits, which are analogous to the 
a bonds (Fig. |5| (b)). In the 1^ band, each valence neutron contains both the positive parity component like a and 
the negative parity one similar to tt. Since the parity of the total system is negative in the 1~ band, the states after 
the parity projection have one neutron in a orbit and the other neutron in tt orbit. Roughly speaking, the 0^, 1~ 
and Oj states are understood as 2a and two valence neutrons in tt^, an and orbits, respectively. The interesting 
point is that the valence neutrons play important roles to develop the clustering structure in the excited bands. We 
can argue that the clustering develops in the 1~ band and mostly in 02" band owing to the cr orbits of the valence 
neutrons, because the a orbit prefers the prolately deformed system as to gain its kinetic energy. This idea originates 
from the application of the two centers shell model to the dimer model for Be isotopes by W. von Oertzen |8| and 
consistent with the argument in the work with the method of AMD+HF by Dote et al. p] . 



C. mechanism of clustering development 

As mentioned above, one of the viewpoints for the mechanism of the clustering development is the molecule-like 
orbits a and tt. Here we try to understand the mechanism from the other viewpoint of the two-center clustering 
model. 

The reason for the clustering development in the ordinary nuclei has been understood as the system gains the 
kinetic energy with the development of clustering although the potential energy looses. In order to understand the 
mechanism of clustering in light unstable nuclei, we think it is helpful to investigate the competition of the kinetic 
and the potential energies in the molecule-like states of ^"Bc. 
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FIG. 4. Density distribution of single-paxticle wave functions of the valence two neutrons in the intrinsic states of Of , 1~ 
and Oj" (left column). The method to extract the single-particle wave functions are explained in the text. The middle and 

right columns arc for the density of the positive parity and negative parity components projected from the single-particle wave 
functions, respectively. The wave functions projected into the parity eigen states are normalized for presentation. 
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FIG. 5. Schematic figures of the molecular orbits tt bond (a) and a bond (b) surrounding 2a clusters. 
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In the density distribution of VAP results, the structures of the 0^, 1~ and states seem to be the two-center 
clustering structures which consist of ^He+a (see figure ^). To estimate the dependence of the kinetic and potential 
energies on the degree of the spatial clustering development we represent the three kinds of configurations for the 
^He+a system corresponding to the 0^, 1~ and states of ^°Be by the simplified AMD wave functions ^amdC^) 
as follows. The intrinsic spin of the single-particle wave functions are fixed to be up or down for simplicity. For 
the ^He+a system with the inter-cluster distance d (fm), the centers of single-particle Gaussian wave functions are 
located around two points ai = {—3d/5y/v,0,0) and 02 = {2d / 5^/^,0,0). The three kinds of configurations of the 
centers for 0^,1^ and states are shown in the right figures of Fig. ^. We put the centers for p f, p |, n t, n i at 
the point di and p "f , p | at 02. The centers for the last 4 neutrons are located at the points very close to a2 as 02 ± (5 
where enough small S is chosen so that the angle 9 between S and 02 is 9— tt/2, n/A and corresponding to the 0^ , 
1~ and 0^ states, respectively. For the parity eigen states projected from these three kinds of the simplified AMD 
wave functions ^amdCZ) we calculate the expectation values of total, kinetic and potential energies as the function 
of the inter-cluster distance d: 



. . ^ ((1 ± P)^AMDm\H\{l ± P)^aMd{'Z.)) 

^ ((1±P)$^md(Z)|(1±P)$amd(Z)) 

, > ^ ((1 ± P)<Pamd{Z)\T\{1 ± P)^amd{Z)) 
^ ((1±P)$amd(Z)|(1±P)$amd(Z)) 

, > ^ ((1 ± P)^amd{Z)\V\{1 ± P)^amd{Z)) 

^ '- ((1±P)$^MI3(Z)|(1±P)$AMI?(Z)) ' 



(21) 
(22) 
(23) 



where we omit the spin-orbit and the Coulomb forces for simplicity. In Fig. ^ we present the total energy, the kinetic 
energy and the potential energy for the three kinds of clustering states 0^, 1~ and 0^ as the function of the distance 
d between ^He and a clusters. As shown in the figure for the total energy (i?), the optimum distances indicate that 
the clustering structure develops in the system for the 1^ state and is most remarkable in the state for Oj, which 
is consistent with the present results of VAP calculations. The shift of the minimum point of the total energy is 
understood by the energetically advantage of the kinetic part as follows. It is found that the kinetic energies in the 
small d region are sensitive to the configurations, while in the case of the potential energies significant differences are 
not seen in the three configurations. When the cluster approaches each other the kinetic energy in the configuration 
for the 1~ state becomes larger than 0^ by l/2huj (about 7 MeV for v = 0.17), because all the neutrons in Qh[ are 
in Os and Op shells, while in the negative parity state for 1~ one valence neutron must rise to the higher sd shell in 
the small d region. Since in the small d limit the wave function is almost same as the harmonic oscillator shell model 
wave function due to the antisymmetrization, the kinetic energy of the sd shell is larger than by l/2hLu than the one 
of p shell. Thus the system for 1~ loses as to the kinetic energy in the small distance d, as a result, the minimum 
point of total energy shift to the larger d region than the one in the 0^ state. In the case of the linear configuration 
for the excited Oj state, the two valence neutrons occupy the sd orbit in the small d limit and the kinetic energy is 
larger than 0^ by Ticu. That is why the optimum point d in 0^ is the largest of the three. In other words, when the 
cluster approaches each other, the system feels the repulsive force in kinetic part because of the Pauli principle. That 
is the reason why the clustering structure remarkably develops in the Oj state. In the analysis with this simplified 
two cluster model, we can conclude that the clustering develops so as to gain the kinetic energy. It is compatible with 
the viewpoint of the molecular a orbit. 



VI. SUMMARY 



We studied the structures of excited states of ^°Be by performing variational calculations after the spin-parity 
projection in the framework of antisymmetrized molecular dynamics (AMD). We explained the formulation of the 
variation after projection and the method to construct the higher excited states with AMD. 

The excitation energies of many levels of ^°Be were reproduced by AMD calculations with the finite range in- 
teractions. The theoretical results well agree to the experimental data of El, E2 transition strength. The results 
successfully describe the new data of B{GT) for the (3 transition strength reduced from the 0° angle cross sections of 
the charge exchange reactions ^''B(t,^He)^°Be*. 

We discuss the structures of excited states. By analyzing the intrinsic structures it is found that the excited levels 
are classified into the rotational bands as K'^ = Of, 2+, 0^ and 1~. Although the a clusters are not assumed in 
the model, the theoretical results show that the 2a + 2n structures appear in most of the states. Particularly the 
molecular-like states with the developed 2a structure constructs the rotational bands = Oj and 1^ in which the 
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FIG. 6. Schematic figures for the intrinsic structure of 0^, 1^ and O2 are shown in the left columns of (a), (b) and (c), 
respectively. Right columns indicate the three kinds of configurations for the centers of Gaussians in the simplified AMD 
wavefunctions which correspond to the two-center ^He+a cluster model for the excited states 0]^(a), l~(b) and Oj(c) of ^"Be. 
The black(white) circles correspond to the centers of Gaussians of the single-particle wave functions for protons(neutrons). 

level spacing is small because of the large deformations. The molecule-like structure develops more largely in K'^ = 1~ 
than in K'^ = Of and most remarkably in = 

We extracted the single-particle wave functions to discuss the behavior of the valence neutrons. In the analysis of 
the intrinsic states of Oj^, 1^ and 0^ states (the band heads of the rotational bands = Oj^, 1^, 0^), we fomid that 
the two valence neutrons are in the molecular orbits surrounding 2a. The positive and negative parts of the orbits 
are analogous to the a and tt bonds in terms of the electron orbits in molecules. Roughly speaking, the Of states, 
1~ and O2 bands is regarded as the states with the two valence neutrons in tt^, aw and cr^ orbits, respectively. The 
neutrons in the a orbit plays as important role for the development of clustering because a orbit profits from the 
kinetic energy in the prolately deformed system. 

In order to understand the mechanism of the clustering development in the excited states of ^°Be, we analyze the 
kinetic and the potential parts of the total energies with the simplified model of the two-center ^He+a cluster model. 
In the 1^ and Of states, when the distance between clusters becomes small the systems feel the repulsive kinetic 
energy because the valence neutrons must rise into higher orbits. Because of the repulsive kinetic energy in the small 
inter-cluster distance the clustering structure develops in the excited states 1~ and Of to gain as to the kinetic energy. 
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